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Abstract 



Diffractive scattering of j*p X + N , where N is either a proton or a nucleonic 
system with Mj\f < 4 GeV has been measured in deep inelastic scattering (DIS) at 
HERA. The cross section was determined by a novel method as a function of the 7*p 
cm. energy W between 60 and 245 GeV and of the mass Mx of the system X up to 
15 GeV at average 

Q2 values of 14 and 31 GeV^. The diffractive cross section da'^^^f /dMx 
is, within errors, found to rise linearly with W . Parameterizing the W dependence by 
the form da'^^-f^/dMx oc (1^2)(2'^ip~2) the DIS data yield for the pomeron trajectory 
cE^ = 1.23 lb 0.02{stat) it 0.04(sysi) averaged over t in the measured kinematic range 
assuming the longitudinal photon contribution to be zero. This value for the pomeron 
trajectory is substantially larger than cF^ extracted from soft interactions. The value of 
cF^ measured in this analysis suggests that a substantial part of the diffractive DIS cross 
section originates from processes which can be described by perturbative QCD. From 
the measured diffractive cross sections the diffractive structure function of the proton 
F^^^\(3,Q'^,Xjp) has been determined, where (3 is the momentum fraction of the struck 
quark in the pomeron. The form F^^^^ = constant ■ (l/xjp)"" gives a good fit to the data 
in all P and intervals with a = 1.46 ± OM{stat) ± OM{syst). 
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1 Introduction 



In deep- inelastic electron-proton scattering (DIS), e~p — > e~ + anything (Fig. |l]), a new class 
of events was observed by ZEUS |l|, 0, ^ and HI characterized by a large rapidity gap 
(LRG) between the direction of the proton beam and the angle of the first significant energy 
deposition in the detector. The properties of these events indicate a diffractive and leading 
twist production mechanism. The observation of jet production demonstrated that there is 
a hard scattering component in virtual-photon proton interactions leading to LRG events. A 
comparison of the energy flow in events with and without a large rapidity gap showed that in 
LRG events the QCD radiative processes are suppressed. 

The diffractive contribution to the proton structure function F2 was measured by HI |^ and 
ZEUS 10. The diffractive electron-proton cross section was found to be consistent with factoris- 
ing into a term describing the flux of a colourless component in the proton and a term which 
describes the cross section for scattering of this colourless object on an electron. LRG events 
were also observed in photoproduction ||^, H . A combined analysis of the diffractive part of the 
proton structure function F2 and the diffractive photoproduction of jets indicated that a large 
fraction of the momentum of the colourless object carried by partons is due to hard gluons 0. 

One of the most interesting questions raised by these LRG events is the precise W dependence 
of the cross section for diffractive scattering of virtual photons on protons, 7*p Xp. Here, W 
is the 7*p cm. energy and the comparison should be done at fixed mass squared of the virtual 
photon, —Q^. In the Regge picture, the elastic and diffractive cross sections in the forward 
direction are expected to behave as (see e.g. [|10[]): 

rfa(t = 0)Moc(iy2)2«^(°)-2 , (1) 

where t is the square of the four-momentum transferred from the virtual photon to the incoming 
proton. From elastic and total cross section measurements for hadron-hadron scattering the 
intercept ajp(O) of the pomeron trajectory ajp{t) was found to be 1.08 |Tl|. A similar energy 



dependence was observed for diffractive dissociation in hadron-hadron scattering of the type 



hih2 hiX, for a fixed mass Mx of the diffractively produced system X (see e.g. [|T^, |T3l)- For 



DIS, with dominantly hard partonic interactions, the BFKL formalism |]T^ leads to a pomeron 



intercept of ajp{0) ^ 1 -|- (12 ln2)as/7r ^ 1.5 at Q'^ = 20 GeV^ which could imply a rapid rise 
of the diffractive cross section with W [T^, T^]. 

In the previous determinations of the diffractive structure function, the subtraction of the 
nondiffr active contribution relied on specific models §[. In the present analysis the separation 
is based on the data. The diffractive contr 

ibution is extracted by a new method which uses the mass Mx of the system X, measured in the 
detector, to separate the diffractive and nondiffractive contributions. The distribution in In 
exhibits, for the nondiffractive component, an exponential fall-off towards small lnM|- values, 
dN"'"'"'^^^^ /din Mx oc exp(6 lnM|^), a property which is predicted by QCD-based models for 
nondiffractive DIS (see e.g. |jT7| , 0). The parameter b of this exponential fall-off is determined 
from the data and is assumed to be valid in the region of overlap between the diffractive and 
nondiffractive components, so allowing subtraction of the nondiffractive background. 

The cross section for diffractive production by virtual photons on protons, 7*p XN, is 
determined integrated over t. The system is either a proton or a nucleonic system with mass 
Mn < 4 GeV. The 4 GeV mass limit results from the acceptance of the detector. 



1 



The prime goal of this analysis is the determination of the W dependence of the diffractive 
7*p cross section in the range 60 < < 245 GeV, Mx < 15 GeV and 10 < < 56 GeV^. 
The paper begins with a brief introduction to the experimental setup and the event selection 
procedure followed by a description of the determination of the mass Mx- Using the measured 
Mx distributions, the widely different behaviour of the nondiffr active and diffractive contribu- 
tions is demonstrated: production of events with low masses Mx is dominated by diffractive 
scattering while nondiffr active events are concentrated at large Mx values. These observations 
lead to a straightforward procedure for extrapolating the nondiffr active background into the 
low mass region and extracting the diffractive contribution. An unfolding procedure is used to 
correct the resulting number of diffractive events in (Mx,VF,Q^) bins for detector acceptance 
and migration effects. From the corrected number of events the cross sections for diffractive 
production by virtual-photon proton scattering are obtained and the W dependence of diffrac- 
tive scattering is determined. Finally, the cross sections are analyzed in terms of the diffractive 
structure function of the proton. 



2 Experimental setup 

The experiment was performed at the electron-proton collider HERA using the ZEUS detector. 
The analysis used data taken in 1993 where electrons oi Eg = 26.7 GeV collided with protons of 
Ep = 820 GeV. HERA is designed to run with 210 bunches in each of the electron and proton 
rings. In 1993, 84 paired bunches were filled for each beam and in addition 10 electron and 6 
proton bunches were left unpaired for background studies. The integrated luminosity was 543 
nb^^. Details on the operation of HERA and the detector can be found in . 



2.1 ZEUS detector 

The analysis relies mainly on the high-resolution depleted-uranium scintillator calorimeter and 
the central tracking detectors. The calorimeter covers 99.7% of the solid angle. It is divided 
into three parts, forward (FCAL) covering the pseudorapidity [| region 4.3 > t] > 1.1, barrel 
(BCAL) covering the central region 1.1 > rj > —0.75 and rear (RCAL) covering the backward 
region —0.75 > t] > —3.8. Holes of 20 x 20 cm^ in the center of FCAL and RCAL are 
required to accommodate the HERA beam pipe. The calorimeter parts are subdivided into 
towers of typically 20 x 20 cm^ transverse dimensions, which in turn are segmented in depth 
into electromagnetic (EMC) and hadronic (HAG) sections. To improve spatial resolution, the 
electromagnetic sections are subdivided transversely into cells of typically 5 x 20 cm^ (10 x 
20 cm^ for the rear calorimeter). Each cell is read out by two photomultiplier tubes, providing 
redundancy and a position measurement within the cell. Under test beam conditions ||20|| , 
the calorimeter has an energy resolution, as, given by oe/E = 1^%/\/E for electrons and 
geIE = 35%/ for hadrons, where E is in units of GeV. In addition, the calorimeter cells 
provide time measurements with a time resolution below 1 ns for energy deposits greater than 
4.5 GeV, a property used in background rejection. The calorimeter noise, dominated by the 



^The ZEUS coordinate system is right-handed with the Z axis pointing in the proton beam direction, here- 
after referred to as forward, and the X axis horizontal, pointing towards the center of HERA. The pseudorapidity 
J] is defined as — ln(tan |), where the polar angle d is taken with respect to the proton beam direction from the 
nominal interaction point. 
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uranium radioactivity, in average is in the range 15-19 MeV for electromagnetic cells and 24- 
30 MeV for hadronic cells. The calorimeter is described in detail in . 



Charged particle detection is performed by two concentric cylindrical drift chambers, the vertex 
detector (VXD) and the central tracking detector (CTD) occupying the space between the 
beam pipe and the superconducting coil of the magnet. The detector was operated with a 
magnetic field of 1.43 T. The CTD consists of 72 cylindrical drift chamber layers organized into 
9 superlayers [^. In events with charged particle tracks, using the combined data from both 



chambers, resolutions of 0.4 cm in Z and 0.1 cm in radial direction in the XY plane are obtained 
for the primary vertex reconstruction. From Gaussian fits to the Z vertex distribution, the rms 
spread is found to be 10.5 cm in agreement with the expectation from the proton bunch length. 

The luminosity is determined by measuring the rate of energetic bremsstrahlung photons pro- 
duced in the process ep — >■ ep7 p2| . The photons are detected in a lead-scintillator calorimeter 



placed at Z = —107 m. The background rate from collisions with the residual gas in the beam 
pipe was subtracted using the unpaired electron and proton bunches. 



2.2 Kinematics 

The basic quantities used for the description of inclusive deep inelastic scattering 

e(A;) +p(P) — > e(A;') + anything 

are: 

= -q^ = -{k-}^f , (2) 



2P-q 
P-q 



(4) 



P-k ' 

= ^^^i^-fM2^^forx«l , (5) 

where k and k' are the four-momenta of the initial and final state electrons, P is the initial 
state proton four-momentum, Mp is the proton mass, y is the fractional energy transfer to the 
proton in its rest frame and W is the 7*p cm. energy. For the range of and W considered 
in this paper we also have W"^ ^ y ■ s, where s = AE^Ep is the square of the ep cm. energy, 
= 296 GeV. 

For the description of the diffractive processes, 

ep e + X + N, 
in addition to the mass Mx, two further variables are introduced: 
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In models where diffraction is described by the exchange of a particle-hke pomeron, is the 
momentum fraction of the pomeron in the proton and (3 is the momentum fraction of the struck 
quark within the pomeron. 



The kinematic variables x, and W were determined with the double angle (DA) method |]23 



in which only the angles of the scattered electron [6'^ and the hadronic system {'~)h) are used. 
This reduces the sensitivity to energy scale uncertainties. The angle characterizes the 
transverse and longitudinal momenta of the hadronic system. In the naive quark-parton model 
7ii/ is the scattering angle of the struck quark. It was determined from the hadronic energy flow 
measured in the calorimeter. A momentum vector ph = {px,Py,Pz) was assigned to each cell h 
with energy E in such a way that p| = -E^. The cell angles were calculated from the geometric 
center of the cell and the vertex position of the event. The angle 'Jh was calculated according 
to 

_ ij:hPxy + {EhPY?-{j:h{E-pz)r . . 

(E.Px)2 + (E.PY)2 + (E.,(i?-Pz)P' ^ ^ 

where the sums, J^h^ here and in the following, run over all calorimeter cells h which were 
not assigned to the scattered electron. The cells were required to have energy deposits above 
60 MeV in the EMC section and 110 MeV in the HAC section and to have energy deposits above 
140 MeV (160 MeV) in the EMC (HAC) sections, if these energy deposits were isolated. The 
last two cuts remove noise caused by the uranium radioactivity which affects the reconstruction 
of the DA variables at low Mx- 

In the double angle method, in order that the hadronic system be well measured, it is necessary 
to require a minimum of hadronic activity in the calorimeter away from the forward direction. 



A suitable quantity for this purpose is the hadronic estimator of the variable y ||2^, deflned by 



We study below events of the type 

ep ^ e + X + rest, 

where X denotes the hadronic system observed in the detector and rest the particle system 
escaping detection through the beam holes. The mass Mx of the system X is determined from 
the energy deposited in the CAL cells according to: 



2.3 Event selection 

The event selection at the trigger level was identical to that used for our F2 analysis ]19[. The 



off-line cuts were very similar to those applied in the double angle analysis of F2 |]T3. The 



resulting event sample is also almost identical to the one used for our recent studies of large 
rapidity gap events in DIS ||^, ^. For ease of reference we list the main kinematic requirements 
imposed, which limit the W and range of the measurement: 
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• i?g > 8 GeV, where E'^ is the energy of the scattered electron, to have rehable electron 
finding and to control the photoproduction background; 

• He < 0.95, where i/e is the variable y calculated from the scattered electron, to reject 
spurious low energy electrons, especially in the forward direction, 

• the impact point of the electron on the face of the RCAL had to lie outside a square 
of side 32 cm centered on the beam axis ("box cut"), to ensure full containment of the 
electron shower, 

• yjB > 0.02, to ensure a good measurement of the angle ■jh and of x, 

• 35 < 5 < 60 GeV, where 6 = J2h (-^ ^Pz), to control radiative corrections and reduce 
photoproduction background. 



The differences with respect to the event selection used for the F2 analysis in |]T9| are an 
increase of the lower limit on E'^ from 5 to 8 GeV and a lowering of the yjB cut from 0.04 to 0.02 
which became possible with the improved noise suppression procedure explained above. The 
increase of the E'^ limit reduces background from photoproduction; the lower yjs cut extends 
the acceptance towards lower W values. It was checked that the F2 values obtained with the 
modified noise procedure were fully compatible with the values published previously in the 
whole Q"^, W range investigated in the present paper. 

The primary event vertex was determined from tracks reconstructed using VXD+CTD infor- 
mation. If no tracking information was present the vertex position was set to the nominal 
interaction point. 

After the selection cuts and the removal of QED Compton scattering events and residual cosmic- 
ray events, the DIS sample contained 46k events. For the analysis of diffractive scattering, 
events with Q"^ > 10 GeV^ were used. The background from beam gas scattering in this sample 
was less than 1% as found from the data taken with unpaired bunches. 



3 Simulation and method of analysis 



3.1 Monte Carlo simulation 



Monte Carlo simulations were used for unfolding the produced event distributions from the 
measured ones, for determining the acceptance and for estimating systematic uncertainties. 

Events from standard DIS processes with first order electroweak corrections were generated 

m for 



with HERACLES 4.4 |2§. It was interfaced using Django 6.0 ^ to ARIADNE 4.03 
modelling the QCD cascade according to the version of the colour dipole model that includes 
the boson-gluon fusion diagram, denoted by CDMBGF. The fragmentation into hadrons was 
performed with the Lund fragmentation scheme |jl8| as implemented in JETSET 7.2 The 
parton densities of the proton were chosen to be the MRSD'- set 



Note that this Monte 



Carlo code does not contain contributions from diffractive 7*p interactions. 
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In order to model the DIS hadronic final states from diffractive interactions where the proton 
does not dissociate, 

— s> e + X + p, 



two Monte Carlo event samples were studied, one of which was generated by POMPYT 1.0 |29 



POMPYT is a Monte Carlo realization of factorizable models for high energy diffractive pro- 



cesses where, within the PYTHIA 5.6 framework, the beam proton emits a pomeron, whose 
constituents take part in a hard scattering process with the virtual-photon. For the quark mo- 
mentum density in the pomeron it has been common to use the so-called Hard POMPYT 
version, (3f{[3) (x [3 ■ {1 — (3). For this analysis the form 

13 f {[3) = constant ■ f3 (11) 

was used which enhances the rate of low Mx events as preferred by the data. 

The second sample was generated following the Nikolaev-Zakharov (NZ) model PT| which was 



interfaced to the Lund fragmentation scheme |32]. In the NZ model, the exchanged virtual 



photon fluctuates into a qq pair or a qqg state which interacts with a colourless two-gluon 
system emitted by the incident proton. In the Monte Carlo implementation of this model the 
mass spectrum contains both components but the qqg states are fragmented into hadrons as 
if they were a qq system with the same mass Mx- Hadronic final states X are generated only 
with masses Mx > 1.7 GeV. For a description of the NZ model see also 0. 

All Monte Carlo events were passed through the standard ZEUS detector and trigger simulations 
and the event reconstruction package. 



3.2 Weighting of diffractive Monte Carlo events 

In order to determine from the number of observed events the number of produced events in 
each bin an unfolding procedure based on a weighted Monte Carlo sample was applied. The 
unfolding procedure is most reliable if the Monte Carlo event distributions are in agreement 
with the data. POMPYT was used for unfolding. However, POMPYT as well as the NZ model 
showed considerable discrepancies relative to the measured W distributions in the kinematic 
range of this study. This problem was overcome as follows: to account for the lack of diffractive 
events in the low mass region, Mx < 1 GeV, events were generated separately for p° production 



via 7*p — i> p°p [|33| and added to the POMPYT event sample. The number of p° events 
and their distribution as a function of W and were determined from the analysis of this 



experiment p3. Furthermore, the POMPYT and p° events were weighted to agree with a 



Triple Regge psj, ^ inspired model (TRM) predicting for the diffractive cross section: 

firj'^^ff ( ]\/f w 

""^m/ = c7.(i + c.-gVMi)(Mo^ + Qr"^(°)(Ml/(Mi + g^)) 

( (12) 

Here C is a normalization constant, afc(O) = (afc(O) = 1) for Mx < Mq {Mx > Mq) and is 
the pomeron trajectory averaged over the square of the four-momentum transfer, t, between the 
incoming and the outgoing proton. The parameters cl, Mq, of the TRM model were deter- 
mined in the unfolding (see section ^ below) and will be referred to as "weighting parameters" . 
The weighted sample of POMPYT events will be referred to as "weighted POMPYT" . 
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3.3 Mass determination 



The mass of the system X was determined from the energy deposits in the calorimeter using 
Eq. |10|. The mass Mjp^"* measured in this way has to be corrected for energy losses in the 
inactive material in front of the calorimeter and for acceptance. The correction was determined 
by comparing for Monte Carlo (MC) generated events the MC measured mass, Mx'""^'^"'^, to 
the generated mass, Mx'"^^"', of the system X. The mass correction was performed in two 
steps. In the first step an overall mass correction factor was determined. In the second step 
the diffractive cross sections were determined by an unfolding procedure (see section |^) taking 
into account for each {Mx, W, Q"^) interval the proper mass correction as determined from the 
MC simulation. 

The overall correction factor f{Mx) was determined from the average ratio of measured to 
generated mass Mx, 

A/fMC'meas 

f{ l\/rMCmeas\ _ ^^^X 

jyiVlx ) — ^MCgen ' 

as a function of Mx,W and Q^. The dependence of f{M^^~) on Mx, W, was found to be 
sufficiently small (±6%) for Mx > 1.5 GeV so that it could be neglected in the first step of the 
mass correction. The average correction factor was /(Mj^^'"™^"*) = 0.68. The same correction 
factor was used for masses below 1.5 GeV. The correction factor / = 0.68 was applied to obtain 
from the measured mass the corrected mass value, M™*" = Mx^°-^ j j ■ 

Figures ^,b show, for MC events, the corrected versus the generated Mx- The error bars in 
Fig. ^3 give the rms resolution for a single Mx measurement. A tight correlation between 
corrected and generated mass is observed except when Mx < 2 GeV where the mass resolution 
is comparable to the value of the mass. The mass resolution increases smoothly from 1 GeV 
near the p mass to 1.3 GeV (2.9 G eV) at M x = 3 GeV (15 GeV). For Mj^^^^"" > 3 GeV it can 

be approximated by a{M^^^'°^) / ^M^^ = 0.75 GeV^/^. 

A test of the MC predictions for the mass measurement at low Mx values was performed by 
studying the reaction 

ep e + p° + p , 

where the pions from the decay p° —>■ tt~^7c~ were measured with the central tracking detec- 



tor |3^. The ir'^TT mass resolution from tracking was 25 MeV(rms). From a total of 60 events 



with 700 < M*o < 800 MeV in the kinematic range Q'^ = 7-25 GeV^, W = 60-134 GeV, 
an average Mx^ of 1.2 GeV and an average mass resolution a{Mx^) of 0.9 GeV were obtained; 
all but 4 events were reconstructed with a mass M^*" below 3.0 GeV. The Monte Carlo sim- 
ulation for this channel predicted M^^^^°^ = 1.2 GeV and a(M|^^^°'^) = 0.8 GeV, in good 
agreement with the data. 

All Mx results presented below refer to Mx^. 



3.4 Acceptance for diffractive events 

A measure of the acceptance for diffractive events is the ratio TZ = j\/^^c:meas ^ j^MCgen 
events measured to events generated in an {Mx, W, Q^) bin using M^^. Figures |^ and ^ show, 
for weighted POMPYT events, the distributions of AT^^^^" (histograms) and (solid 
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points) for the (W, Q"^) bins used in this analysis (see section . Here, the generated values 
for MxiW^Q"^ were used for J\f^'^(^9<^^ while and the double-angle quantities for W and 

were used for j\fMCmeas^ ^.j^g analysis of the data. The Mx distributions increase from 
small Mx values to a maximum at Mx = 2 — 5 GeV and then fall off towards higher masses. 
There is some leakage of events into the low Mx bin as seen in the ratio TZ shown in the 
right-hand parts of Figs. ||, ^. 

The shaded areas mark the Mx regions used for extraction of the diffractive cross sections. The 
ratio TZ is above unity at small Mx as a result of the migration from higher Mx masses; for 
larger Mx values TZ is rather constant and between 70 and 100% in the bins considered for the 
analysis, except in the highest W interval for = 14 GeV^ where the acceptance is around 80% 
at low masses falling to about 50% at Mx = 15 GeV. This is caused by the reduced efficiency 
for detecting the scattered electron and by the requirement that 6 = J2h{E — pz) > 35 GeV. 

3.5 General characteristics of the Mx distributions 

The method of separating the diffractive and nondiffr active contributions is based on their very 
different Mx distributions. As a first illustration. Fig. ^ shows the distribution of Mx versus 
W for the data. Two distinct classes of events are observed, one concentrated at small Mx, the 
second extending to large values of Mx- Most of the events in the low Mx region exhibit a large 
rapidity gap, which is characteristic of diffractive production. This is shown by Fig. |^ where 
the events with a large (small) rapidity gap, r]max < 1-5 ( rjmax > 1-5) are marked by different 
symbols. Here rjmax is the pseudorapidity of the most forward going particle. For this analysis 
a particle is defined as an isolated set of adjacent calorimeter cells with more than 400 MeV 
summed energy, or a track observed in the central track detector with more than 400 MeV 
momentum. A cut of rimax < 1-5 corresponds to a visible rapidity gap larger than 2.2 units 
since no particles were observed between the forward edge of the calorimeter [rj = 3.7 — 4.3) 
and rj = 1.5. For rjmax < 1-5 the contribution from nondiffractive scattering is expected to be 
negligible |], §. 

The measured Mx distributions are shown in Figs. P(a-c) for three W intervals, = 90 — 
110, 134 - 164 and 200 - 245 GeV at = 14 GeV^. The distributions are not corrected for 
acceptance. For all W bins two distinct groups of events are observed, one peaking at low Mx 
values, the other at high Mx values. While the position of the low mass peak is independent 
of W, the high mass peak moves to higher values as W increases. As already seen, most events 
in the low mass peak possess a large rapidity gap. This is illustrated by the shaded histograms 
which represent the events with rjmax < 1-5. 

The size of the rapidity gap, Ar/, can be seen from Figs, ^(g - i) which show the distributions 
of the rapidity gap between the edge of the calorimeter {rj ^ 3.9, which is the rj value of the 
geometric center of the HAG cells closest to the proton beam) and the most forward lying cell 
with energy deposition greater than 200 MeV (the threshold is reduced in comparison with the 
determination of rjmax because here single cells are considered instead of a group of cells). The 
plots give the distribution of all events (points with error bars) and those with Mx < 3 GeV 
(shaded) and 3 - 7.5 GeV (skewed hatching). There is a strong concentration of events at small 
rapidities, Arj < 1, which stem from non diffractive processes. The distributions demonstrate 
that the majority of low Mx events are associated with a large rapidity gap Arj > 2. The 
average rapidity gap increases with growing W. For a given W value rjmax is correlated with 
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the maximum possible Mx value but does not allow the determination of Mx uniquely because 
of fragmentation effects. Since our aim is the determination of the diffractive cross section as 
a function of Mx, the analysis was based on Mx and not on rjmax- 

In Fig. 1^ the measured Mx distributions are compared with the NZ and CDMBGF predictions. 
The shaded distributions show the NZ predictions for diffractive production. They peak at 
small masses. The predictions of CDMBGF for nondiffractive production (dotted histograms) 
peak at high masses. The sum of the NZ and CDMBGF contributions reproduce the main 
features of the data which are the low and high mass peaks. 

The properties of the Mx distributions can be understood best when plotted as a function of 
InM^, shown in Figs, ^(d-f) and Figs. |^(d-f). Here, and in the following, masses and energies 
are given in units of GeV. In this representation the low mass peak shows up as a plateau- 
like structure at low lnM|^, most notably at high W values. The high mass peak exhibits a 
steep exponential fall-off towards smaller In Mx values. The shape of the exponential fall-off 
is independent of W , a property which is best seen when the XnM'j^ distributions are replotted 
in terms of the scaled variable [lnM|- + ln(s/H^^)] (the total ep cm. energy squared, s, is 
introduced for convenience). This is shown in Fig. ^ at = 14 and 31 GeV^ where the scaled 
In Mx distributions are overlaid for three W intervals. The position of the high mass peak in 
In Mx grows proportionally to In and the slope of the exponential fall-off towards small 
In Mx values is approximately independent of W . 



3.6 Mx dependence of the nondiffractive contribution 

While in diffractive scattering the outgoing nucleonic system remains colourless, in nondiffrac- 
tive DIS the incident proton is broken up and the remnant of the proton is a coloured object. 
This gives rise to a substantial amount of initial and final state QCD radiation, followed by 
fragmentation, between the directions of the incident proton and the current jet as illustrated 
in Fig. |l[ The salient features of the resulting Mx distribution, namely the exponential fall-off 
and the scaling in [lnM|^ -|- ln(s/Vr^)], can be understood from the assumption of uniform, 



uncorrelated particle emission in rapidity y along the beam axis in the 7*j9 system |37 



dNpart/dy = \, A = constant. (13) 

At the values under study, for 3^ > 0, the beam axis in the 7*p system is approximately 
given by the proton direction in the HERA system. Since the shape of the rapidity distribution 
is invariant under translations along the 7* beam axis we translate the y distribution measured 
in the 7*p system until the point of maximum rapidity agrees with the maximum rapidity ymax 
in the ep system. For an (idealized) uniform y distribution between maximum and minimum 
rapidities of ymax and 3^mm, respectively, the total center of mass energy W is given by 

W"^ = Co ■ exp{ymax - ymin), aSSUmiug (3^max - 3^mm) > 1- (14) 

Here, cq is a constant. The mass Mx of the particle system that can be observed in the detector 
is reduced by particle loss mainly through the forward beam hole: 

Mi = Co ■ exp(3;,t™* - ymin) = W ■ expiyf:^,, - ymax) (15) 

where yf^it denotes the hmit of the FCAL acceptance (neglecting the mass and transverse 
momentum of the produced particles). Equation ^ predicts scaling of the InM^ distribution 
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when plotted function of \n{M'^/W'^), in agreement with the behaviour of the data seen in 
Fig. H where the data are plotted in terms of [lnM|. + ln(s/Vr^)]. The quantity {ymax — yumit) 
the effective width of the beam hole and can be estimated from the effective maximum rapidity 
and the detector geometry. 

The mass Mx is expected to fluctuate statistically due to a finite probability -P(O) that no 
particles are emitted between yu^it yf^n — A3^. This generates a gap of size A3^. The 
assumption of uncorrelated particle emission leads to Poisson statistics which predicts -P(O) = 
exp{—XAy) resulting in an exponential fall-off of the InM^ distribution, 

^j\J'nondif f 

^^^=cexp(M„Mj) , (16) 

where the slope b is equal to the parameter A and c is a constant. The exponential fall-off of the 
In Mx distribution towards small values of In Mx expected from this simplified consideration is 
indeed observed in models which include QCD leading order matrix elements, parton showers 
and fragmentation such as CDMBGF shown by the dashed histograms in Figs. |^(d-f). This 
shows that the known sources of long range correlations like conservation of energy and momen- 
tum, of charge and of colour, which are incorporated in CDMBGF, do not lead to significant 
deviations from an exponential behaviour with the possible exception of a very small fraction 
(0.2 — 0.4%) of the CDMBGF events which is found above the exponential at low lnM|^ values. 

In principle, the exponential fall-off of the InMx distribution should start at the maximum 
value of InM^ allowed by kinematics and acceptance, Max(lnM^) = InW^^ — (2 to 3). The 
data in Fig. |^d-f and Fig. ^ break away from the exponential behaviour towards high values of 
In Mx leading to a rounding-ofF. It mainly results from the finite size of the selected W intervals, 
the edge of the calorimeter acceptance in the forward direction (r/edge = 3.7 to 4.3) and the finite 
resolutions with which W and Mx are measured. With good accuracy the exponential fall-off 
is observed for lnM|^ < InVT^ — rjo, with r]o ~ 3.0, over more than two units of rapidity (see 
also Section 4.2). 

We would like to add three remarks. Firstly, the value of the slope b is little affected by detector 
effects: it is almost the same at the detector level as at the generator level. This was verified by 
Monte Carlo simulation of nondiffr active events with CDMBGF. The MC events were selected 
using the same selection cuts as for the data. The mass Mx of a standard nondiffractive DIS 
event at the generator level was defined as the invariant mass of all particles (excluding the 
scattered electron) generated with pseudorapidities rj < 4.3, the nominal end of the detector. 
In Fig. 1^ the MC generated (dashed histograms) and MC measured mass distributions (solid 
histograms) are shown. The exponential slope values of6 = 1.9±0.1 (shown as straight lines) 
obtained at the generator level agree with the slope values found at the detector level to within 
±0.1 units compatible with the statistical errors of the simulation. 

Secondly, the value predicted by CDMBGF for the slope is 6 = 1.9 ± 0.1 while the data 
yield a shallower slope, the average being b = 1.46 ± 0.15. Note, the value of the slope b 
cannot be predicted precisely by the models; rather, DIS data have to be used to fix the 
relevant parameters of the CDMBGF model. For instance, short range correlations arising from 
resonance production affect the value of b. The observed difference between the exponential 
fall-off found in the data and predicted by the Monte Carlo simulation indicates that estimating 
the tail of the nondiffractive background from Monte Carlo simulation alone may lead to an 
incorrect result for the diffractive cross section. 
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Thirdly, we determine the slope b from the data in a region where the nondiffr active contribution 
dominates. The exponential fall-off will be assumed to continue with the same slope into the 
region of overlap with the diffractive contribution. The nondiffractive event sample generated 
with CDMBGF indicates a small excess of events above the exponential fall-off (see above). If 
allowance is made for a similar deviation from the exponential fall-off in the data, the numbers 
of diffractive events obtained after subtraction of the nondiffractive background change by less 
than ^ 30% of their statistical error. 



3.7 Mx dependence of the diffractive contribution 

In diffractive events, the system X resulting from the dissociation of the photon is, in general, 
almost fully contained in the detector while the outgoing proton or low mass nucleonic system, 
j^dtssoc^ escapes through the forward beam hole. Furthermore, diffractive dissociation prefers 
small Mx values and leads to an event distribution of the form dM /dMx oc 1/(M|^)^^+"^ or 

' , (17) 



din Ml (m; 



approximately independent of W . At high energies and for large Mx one expects n leading 
to a roughly constant distribution in lnM|^. Such a mass dependence is seen in diffractive 



dissociation of pp scattering (see e.g. ||T2|, |T3|)- A value of n ~ is also expected in diffractive 



models as the limiting value for the fall-off of the mass distribution (see e.g. the NZ model 



To summarize this and the previous section, the diffractive contribution is identified as the 
excess of events at small Mx above the exponential fall-off of the nondiffractive contribution 
with lnM|^. The exponential fall-off permits the subtraction of the nondiffractive contribution 
and therefore the extraction of the diffractive contribution without assuming the precise Mx 
dependence of the latter. The Mx distribution is expected to be of the form 

D + c exp(6 In M|), for In M^ < In VT^ - r/o. (18) 



d \nMl 



Here, D denotes the diffractive contribution, the second term represents the nondiffractive 
contribution and InW^^ — r/o is the maximum value of lnM|^ up to which the exponential 
behaviour of the nondiffractive part holds. We shall apply Eq. |T^ in a limited range of In M^ 
for fitting the parameters 6, c of the nondiffractive contribution. The diffractive contribution 
will not be taken from the fit result for D but will be determined by subtracting from the 
observed number of events the nondiffractive contribution found in the fits. 



3.8 Contribution from nucleon dissociation 

The contribution from the diffractive process where the nucleon dissociates, 

ep^e + X + N'^i''"^, 

was estimated by Monte Carlo simulation. Assuming factorisation and a Triple Regge formalism 
351 p6| for modelling 7*p — > X+p and 7*]? — > X-|-iV'^*'^*°'^, the measured cross sections for elastic 
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and single diffractive dissociation in pp (and pp) scattering, pp p + p and pp ^ p + ]\[dissoc^ 
were used to relate 7*p ^ X + ]S[dissoc ^.^ x + p. 

The secondary particles from jsj^issoc ^g^g^y -^gj^g found to be strongly coUimated around the 
direction of the proton beam. Analysis of the angular distribution of the secondary particles 
as a function of showed that for < 2 GeV basically no energy is deposited in the 
calorimeter while for events with Mjy > 6 GeV there are almost always secondaries which 
deposit energy in the calorimeter. Furthermore, events of the type ep —>■ e + X + ]\fdissoc^ 
where decay particles from ]\fdissoc (jgposit energy in the calorimeter, have in general a mass 
reconstructed from all energy deposits in the calorimeter (including those from the decay of the 
system X but excluding those of the scattered electron) which is much larger than the mass 
of X. As a result these events make only a small contribution to the event sample selected 
below for diffractive production of 7*p ^ X + N either with Mx < 7.5 GeV or Mx < 15 GeV. 
To a good approximation, the selection includes all events from dissociation of the nucleon, 
7*p — > X + JV'^**'^°'^, with Mtv < 4 GeV. From the comparison with the pp data, we estimate the 
contribution of 7*p X + A^'^«^'^°'= with < 4 GeV to the diffractive sample to be 11 ± 5%. 



4 Extraction of the diffractive contribution 
4.1 Binning in and W 

The cross section for 7*p —>■ X + N was determined for two intervals, 10 - 20 GeV^ and 
20 - 56 GeV^, the average values being 14 GeV^ and 31 GeV^, respectively. This choice of 
intervals was motivated by the available event statistics and by the requirement of good 
acceptance. These requirements also determined the W range of 60 < < 245 GeV. The 
intervals in W were chosen so as to have equidistant bins in In providing approximately equal 
numbers of events in each W bin. For the bin width, A In = 0.4 was used, commensurate 
with the resolution for Iniy^ of 0.17 for diffractive events and 0.32 for nondiffractive events, as 
determined by Monte Carlo simulation. The better resolution for diffractive events is due to 
the fact that here the decay particles from the system X are almost completely contained in 
the detector. 

The total number of accepted events in the region 60 < W < 245 GeV and 10 < < 20 GeV^ 
(20 < g2 < 56GeV2) was 14466 (11247). 



4.2 Fitting the nondiffractive contribution 

The mass distributions for all W and intervals are presented in Figs. |TD| - |rT| in terms of 
lnM|-. The diffractive contribution was obtained for each W,Q^ interval by determining first 
the nondiffractive background from a fit to the data at large InM]^. The nondiffractive back- 
ground was then extrapolated to the small In Mx region and subtracted from the data giving 



the diffractive contribution as illustrated in Fig. 12. The fit for the nondiffractive background 



was performed using Eq. |T8|. The diffractive contribution D was assumed to be of the form 

A.^. (19) 
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where Di and D2 are constants. The second term allows for contributions from lower lying 
Regge poles contributing to 7* pomeron scattering at a cm. energy of Mx or from a hard 
quark distribution in the pomeron (e.g. Hard POMPYT). The fit parameters are Di, D2, b 
and c. 

The fits were performed to the data in the range 

lnQ2 < lnM|. < Max(lnM|^). The lower 
limit for lnM|^ was chosen according to the expectation of the diffraction models 
that for Mx > the diffractive contribution is of the form given by Eq. |19|. The upper 
limit Max(lnM|.) was chosen as the maximum value of InM]^ up to which the data exhibit 
an exponential behaviour. The maximum value of In Mx was determined by fitting the In Mx 
distributions for each (H^, Q^) interval with a varying maximum value of InM]^. In most 
{W,Q^) intervals a clear boundary as a function of Max(lnM|-) was observed beyond which 
the probability for the fit dropped rapidly. The boundary marks the location where the 
distribution starts to deviate from an exponential behaviour. The boundary was found to be at 
In Mx + ItlIs/W"^) = 8.8 corresponding to a value of 770 = 3.0, a value which is in good agreement 
with expectations of ?7o = 3^m.ax — 3^/1mji (see discussion above). The fits were performed with 
Max(lnM|^) = ln{W'^/s) + 8.6. In the studies of systematic uncertainties the maximum value 
was increased (decreased) by 0.2 (0.4), see section The probabilities for the fits were on 
average 40%. 

The fits were performed by including both terms, Di and D2, (extended fits) as well as by 
assuming D2 = (nominal fits). Since the two fit procedures resulted in only minor differences 
we used the results from the nominal fits and used those from the extended fits for estimating 
the systematic errors (see section ^.2|) . The solid lines in Figs. [1^ - |ll] show for all and 



W bins the exponential fall-off of the nondiffr active contribution resulting from the fits. The 
nondiffr active contribution moves to larger In M|- values proportional to In W"^. As W increases, 
the diffractive contribution appears with little background over an increasing Mx range. Above 
= 90 (164) GeV the nondiffractive background is small up to Mx values of 7.5 (15) GeV 
(see also Table 1 below). 

In Figs. and |rT| the nondiffractive background estimates are also compared with the predic- 
tions of the CDMBGF simulation. The dotted histograms display the predictions of CDMBGF 
for the nondiffractive contribution normalized to 85% of the number of events observed in the 
data in each (W, Q^) bin, the 15% reduction accounting roughly for the diffractive contribution. 
The qualitative features of the data are reproduced by CDMBGF in all and W intervals, 
although the exponential slopes are steeper than in the data (see discussion in section pT^ ). 



5 Determination of the diffractive cross section 

The number of diffractive events, J^mlls^ determined in all and W bins for the Mx 
intervals < 3 GeV (lnM|. < 2.2), 3 - 7.5 GeV (lnM|. = 2.2 - 4.0) and 7.5 - 15 GeV (lnM|. = 
4.0 — 5.4) by subtracting from the observed number of events, Afobs, the contribution from 
electron beam gas scattering, J\f'^9as ^ ^^^^ ^j^g nondiffractive contribution, J\f^°^'^^Sf ^ obtained 
from the fit, M^H^ = Mobs — M'^^"'^ — J\f^°^'^^ff . Electron-gas scattering produces low Mx 
events, which are candidates for diffractive scattering; most of them have no reconstructed 
vertex. Events from unpaired electron bunches were used to obtain A/'^^"''. 
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From M^^as number of produced diffractive events, MplH-, was obtained by a Bayesian 
unfolding procedure which took into account detector effects such as bin-to-bin migration, 
trigger biases and event selection cuts. In the unfolding, the event distribution generated by 
Monte Carlo, ndi), was reweighted as discussed in section so that the observed Monte Carlo 
distribution, no{j), reproduced closely AA^ea^, the event distribution measured in the data which 
will be denoted by rioat^j). The indices i, j denote the three - dimensional bins in (Mx, W , Q"^) 
in which the data were analyzed. For every set of weights we determined from the generated 
distribution the observed distribution and the transfer matrix TGo{j,i), which leads from the 
observed to the generated distribution, ndi) = J2j TgoU, '^)'^o(j)- In an iterative procedure the 
obtained from the differences of nuatii) and no{j) was minimized by varying the weighting 
parameters of the TRM function given in Eq. 12. A good description of the data was obtained; 
the minimum = 34 for 28 degrees of freedom. The values found for the weighting parameters 
were cl = 0.1, Mq/Q = 1 and = 1.2. The resulting matrix TcoUii) was then used to 
determine the unfolded event distribution from that measured, riuii) = J2jTGo{j,i)nDat{,j)- 
The number of unfolded events is denoted by MplH- In the calculation of the statistical errors, 
the bin-to-bin correlations were neglected. The errors were checked by examining the spread 
of results obtained by dividing the data into several subsamples. 



5.1 Evaluation of the cross sections 

For the final analysis, only bins where the fraction of nondiffractive background was less than 
50 % and the purity was above 30 % were kept. Purity is defined as the ratio of the number of 
events generated in the bin and observed in the same bin divided by the total number of events 
observed in the bin. The average purity in a (Mx, VF, Q"^) bin was 43 %. 

Electromagnetic radiative effects were corrected for in every (Mx, Q^) bin [39|. The correc- 
tions were less than 10% and independent of W . 

The average differential cross section for ep scattering, in a {Mx^W^Q"^) bin, is obtained by 
dividing the number of unfolded events, J^prld, by the luminosity and the bin widths. The 
lower limit of Mx, was taken to be 2m^, where is the pion mass. 

The cross sections for the process ep —>■ eXN can be expressed in terms of the transverse (T) 
and longitudinal (L) cross sections cr^"^^ , cr'j^^^ , for 7*p — *■ XN as follows pOf: 



^ dMxdlnW'dQ' ~27r^^^ + IJi^T + J[1 (i _ + i ^..// ^ ^..//J-(^O) 

The relative contribution of the correction term in the third square bracket on the r.h.s. of 
Eq. [2^ is negligible if |/ -C 1 or a'j^^^ <C a^-^^ . Since y ~ W^/s, the contribution can be 
substantial only at high W values. In the extreme case that af'-^^ ^ o'^^^ , the correction 
term will increase [a^-^^ + cr'j^^^] by at most 35% for the highest W bin (200 — 245 GeV) and 
11% for the next highest W bin (164 - 200 GeV). If af^^ = Q'^/Ml ■ a^^^, as e.g. in the 
Vector Dominance Model or in partonic models (see e.g. [0), then the second term increases 
[a^-^^ + o'^j^^^] by at most 31% (17%, 6%) for the bins with the highest W , values and 
Mx < 3 GeV (3 - 7.5 GeV, 7.5 - 15 GeV) faUing below 10% (6%, 2%) in the next highest W 
binQ. 

^Diffractive p° production via ep — > ep°p is a known contribution to cr^*^''^ (see |0). The fraction of 
diffractive events in the lowest Mx bin (2m7r < Mx < 3 GeV) from p° production is around 20 % for the W ^ 
region under study. 
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In the following analysis, the correction term was set equal to one: 



(21) 




a {l-yf + l dMxdlnW^dQ^ 



The numbers of events observed and estimated to come from background are given in Table 1, 
together with the values of the ■y*p — > XN differential cross sections averaged over the specified 
Mx range. The cross sections are quoted for the values of 14 and 31 GeV^ and the W values 
corresponding to the logarithmic means of the W interval hmits. The average Mx values in 
each Mx interval are 1.9, 5.1, 11.0 GeV at = 14 GeV^ and 2.0, 5.1 and 11.0 GeV at = 31 
GeV2. 

5.2 Systematic errors 

The systematic errors on the cross sections were estimated by varying the cuts and algorithms 
used to select the events. The bin-by-bin changes from the standard values were recorded. For 
reference each test is numbered; the number is given in brackets {}. 

The efficiency for finding the scattered electron was around 100% for E'^ > 18 GeV, falling to 
55% at E'^ — 10 GeV. To evaluate the resulting uncertainty on the cross section the cut on E'^ 
was raised {1} (lowered {2}) from 8 to 10 GeV (7 GeV); the box cut was varied from 32 cm to 
36 cm {3} (28 cm {4}). The changes on the measured cross sections were negligible except for 
one low W bin where a 28% change of the cross section was observed. 

In order to test for remaining background from photoproduction and for the sensitivity to 
radiative effects, which were not simulated in the diffractivc Monte Carlo, the cut on J^hi^^Pz) 
was raised {5} (lowered {6}) from 35 GeV to 38 GeV (32 GeV). This resulted in small changes 
of the cross sections, except in one bin where it reached 20%. 

The effect of the cut on yjB was tested by raising {7} (lowering {8}) it from 0.02 to 0.03 (0.01). 
The changes were found to be below 5%, except for the low Mx interval where in two low W 
bins the higher yjB cut substantially reduced the acceptance; for these bins changes of 19% 
and 47%, respectively, were observed. 

The mass correction factor was assumed to be different by -|-10% {9} (-10% {10}) in the Monte 
Carlo simulation and in the data. This affected the low W bins where changes up to 17% were 
observed. 

The analysis was performed without requiring a reconstructed event vertex and could therefore 
be affected by background from beam-gas scattering. The analysis was repeated requiring an 
event vertex as determined by tracking where the Z-coordinate of the vertex had to lie in the 
interval -50 cm to +40 cm {11}. This requirement was satisfied by 89% of the accepted events. 
The vertex requirement reduced the total number of events in the accepted {Mx-,W,Q'^) bins 
obtained with unpaired electron bunches from 12 to 1 event thereby suppressing the electron- 
gas background to a negligible amount. Apart from this effect reductions of the cross sections 
by at most 12(16)% in the highest W bins at = 14(31) GeV^ were observed. For the other 
bins the changes were smaller. 

Uncertainties in the number of diffractive events resulting from the subtraction of the non- 
diffractive background were estimated by increasing {12} (decreasing {13}) the upper limit of 
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the fit range in InAf^ by 0.2 (0.4) units; by decreasing the lower limit of the fit range in lnM|- 
by 0.4 units {14}; by repeating the fits for the nondiffractive background with the extended 
form Eq. |1^, D2 7^ 0, for the diffractive part {15}. The typical changes were below 10%. The 
largest change was observed for the lowest W bin at = 31 GeV^ where it reached 37%. 

The uncertainty resulting from the Monte Carlo weighting procedure was estimated by varying 
the weighting parameter from 1.2 to 1.1 {16} and 1.3 {17}. To estimate the uncertainty of 
the Monte Carlo modelling for the low Mx region the p° events for Mx < 1 GeV were replaced 
by a system decaying into 7r+, tt^, 7r°'s with a phase space like distribution. Changes of less 
than 8% were observed. 

The total systematic error for each bin was determined by adding quadratically the individual 
systematic uncertainties, separately for the positive and negative contributions. The total errors 
were obtained by adding the statistical and systematic errors in quadrature. The errors do not 
include an overall normalization uncertainty of 3.5% of which 3.3% is from the luminosity 
determination and 1% from the uncertainty in the trigger efficiency. 



6 Differential cross section for 'y*p ^ X + N 

The differential cross section da'^JJ^xN/^Mx was determined according to Eq.(^) for the 



different {Mx^V iQ'^) intervals. The results are shown in Fig. ^ as a function of W averaged 
over the Mx bins 2m^ - 3 GeV, 3 - 7.5 GeV and 7.5 - 15 GeV at = 14 and 31 GeV^. The 
solid points show the measured values. The inner error bars give the statistical errors. For the 
full bars the statistical and systematic errors have been added in quadrature. 

The cross section, within errors, is seen to rise linearly with W at both values for all Mx 
bins up to 7.5 GeV. For the Mx bins (7.5 - 15) GeV the data in the accepted W range are 
consistent with this rise. 



In a Regge - type description the W dependence of the diffractive cross section is of 

the form 

d(rJpLxN{Mx,W,Q'^,t) ^ .^2^2a^{0)-2 . gi(B+2a^'ln{TyV(M2,+Q2)) C221 

dtdMx ^ \ ) 

where oijpit) = ajp(O) + a^'t is the pomeron trajectory and B and a^' are parameters. The 
cross sections in each {Mx, Q^) interval were fitted to the form 



dM 



X 



where IT^ stands for ajp(t) averaged over the t distribution. The fit was performed by con- 
sidering cF^ and the six normalization constants for the six {Mx, Q"^) intervals as free param- 
eters. Taking into account only the statistical errors, cF^ was found to be 1.23 ± 0.02 with 
X^/ndf = 11.7/24. The systematic uncertainties were estimated by repeating the fit indepen- 
dently for every source of systematic error discussed above. The results are shown in Fig. |l^. 
The lowest value obtained was 1.20, the highest value was 1.25. The observed deviations 
were added in quadrature leading to the final value: 
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= 1.23 ± 0.02{stat) ±OM{sy St). 

The W dependence was also determined by restricting the analysis to those {Mx ,W jQ"^) bins 
where almost all events have a large rapidity gap and where, therefore, nondiffractive contribu- 
tions should be negligible. Of the events observed with Mx < 3 GeV, 98% have a rapidity gap 
Arj > 2. The fit to these events yielded cF^ = 1.24 ± 0.03 l^os- Note that the average rapidity 
gap for these events grows with W and a possible background from nondiffractive scattering 
would diminish with rising W. Consequently, if there were a nondiffractive contribution left 
in the sample, the correction for this background would lead to an increase of the value of 
compared to the result obtained. 

The effect of the kinematically allowed minimum value of \t\ (|t|mm) on the t - integrated cross 
section and therefore on the value of (T^ is negligible, |t|mm being less than 10^^ GeV^. If 
a^p' = then is equal to the pomeron intercept at t = 0, ajp{0) = cT^. If the slope a J in 
diffractive DIS is the same as for the soft pomeron (a^' = 0.25 GeV~^) and the parameter B is 
equal to half the value observed for elastic pp scattering |jl2| {B = 4.5 GeV~^), then ajp{0) will 



increase from 1.23 to 1.26. Our result can be compared with the soft pomeron trajectory |TT 
C(jp{t) = 1.08 + 0.25t, as determined from hadron - hadron scattering. Assuming B = 4.5 GeV~^ 
and averaging over t, the soft pomeron predicts = 1.05. In extracting the diffractive cross 
section, the assumption a'f^-^^ = was made. Assuming instead a'^^^ = {Q'^/Mx)a^^^ (see 



for example M^) will increase from 1.23 to 1.28. Hence, a positive slope of the pomeron 



trajectory and/or a finite a'^^^ contribution will lead to a larger Q!jp(0) value and increase the 
difference with respect to the soft pomeron intercept. 

The observation of being substantially larger in diffractive DIS than expected for the soft 
pomeron is in line with the expectation of perturbative QCD and shows that deep inelastic 
diffractive scattering has a perturbative contribution. The measured value of is smaller than 
the value which would follow from the BFKL formalism, ajp(O) ~ 1.5. It is in broad agreement 



with the effective pomeron intercept expected in the perturbative models of |^2|, |4^ where the 
Bjorken-x dependence of the gluon momentum density xg{x, Q"^) of the proton determines the 
W dependence of diffractive scattering. 

There are also models where the effective ajp(O) is expected to be smaller in diffractive hadron - 
hadron or photon - hadron scattering as compared to deep inelastic diffractive scattering because 
in hadron - hadron or photoproduction processes, in addition to single, multiple pomeron 
exchanges also contribute (see e.g. [Q). In these models the importance of multiple pomeron 



exchanges decreases quickly with growing photon virtuality. 



6.1 Diffractive contribution to total deep inelastic scattering 

The relative contribution of diffractive scattering to the total virtual-photon proton cross sec- 
tion, was determined for the W bin (164 - 200 GeV) where data on the diffractive cross 
section are available up to Mx = 15 GeV. The data for a*** were taken from the analysis of the 



proton structure function F2 [0. The ratio r'^*-^-^ = / dMxcy^l^J^xNi^N < 4 GeV)/(j**p is given 
in Table 2 integrated over different Mx bins at = 14 and 31 GeV^. In the lowest Mx bin 
the relative contribution from diffractive scattering to the total DIS drops by a factor of about 
three from = 14 to 31 GeV^. With increasing Mx the relative contributions from diffractive 



17 



scattering tend to become equal for the two values of . The observed behaviour does 
not preclude a leading twist behaviour of the diffractive DIS cross section (observed by [|l], H): 
the measurements for the two different values correspond to different values of x, namely 
X = Q'^/W'^ = 4 • 10~^ and 9 ■ 10"'^, respectively. Furthermore, it is conceivable that for fixed 
X the behaviour of the diffractive cross section changes with Mx and only its integral over 
Mx is of leading twist [ |H . 



7 Diffractive structure function of the proton 

The DIS diffractive cross section, ep ^ e + X + N, can be expressed in terms of the diffractive 
structure function F^^^\l3,Q'^,Xjp) as follows p6| : 

da'^;^SeXNif^,Q^XjP,MN <^GeV) _ 2na\ ^ n{3) 



[i + {i-yY]Fr\P.Q\x,) (24) 



if the contribution from longitudinal photons is neglected. 

In Fig. ^ we show F2^^^ of this analysis as calculated from the differential cross sections 
da'!^l^J_^XN/dMx, Mx < 4 GeV. The result is plotted as a function of x^ for different values 
of j3 and (solid points). The error bars show the statistical and systematic errors added 
in quadrature. The data from our previous analysis of F2^''^'* are also shown as the open 
points. In the previous analysis, the diffractive contribution was determined with a rapidity 
gap method using CDMBGF to estimate the nondiffr active part. Note that these data have 
been evaluated at somewhat different (3 and values than in the present analysis. 

The diffractive structure function falls rapidly with increasing x^, the x^ dependence be- 
ing the same within errors in all (3 and intervals. A good fit to the data from this 
analysis (solid points) is obtained with the form ^2^''^'' = constant ■ (l/XpY yielding a = 
1.46 ± 0.04{stat) ± 0.08{syst). Note, for fixed and (3 the x^ dependence of F<^^^^ is equiva- 
lent to the W dependence of da'^^^^ {Yv ~^ XN)/dMx discussed above, the values of a and 
being connected by the relation IT^ = (a + l)/2. 

The value of a measured in the present analysis is somewhat higher than the value of 1.30 ± 
O.OS{stat)toAl{syst) found in our previous F2 analysis which can be understood by the 
way the nondiffractive background was subtracted. To investigate the effect of the different 
background estimates we subtracted the nondiffractive contribution using CDMBGF as in the 
previous analysis and determined a for the same kinematic region in x^ and (3. The result 
was a = 1.28 ± 0.04 while the new method gave a = 1.42 ± 0.08. In the Mx > 3 GeV 
region the difference to the previous analysis is due to the new method of estimating the 
nondiffractive background. For Mx < 3 GeV both methods in this analysis gave the same result 
(a = 1.48 ± 0.06) which is not surprising since in this region the nondiffractive contribution is 
found to be negligible in both methods. 



Figure ^ shows also the ^2^*^^'' values obtained by the HI collaboration |^. HI found 
1.19 ± 0.06 ± 0.07, a value which is smaller than the result obtained in this analysis. 



a 



^The previous analysis was limited to the region Mx > 2.8 GeV. 
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The dependence of F^^^'' on (3 in this analysis was determined as follows: the largest range in 
/? is covered for = 0.003 as can be seen in Fig. |T3|. We chose, therefore, for every (Mx, 
Q"^) bin that F^^^^^ value with closest to 0.003 and determined by interpolation with the 
expression F^^^'^ = (x^/0.003)'' ■ F2^^\x^) the value of F^^^\Q'^, P,x^ = 0.003). The result 
is shown in Fig. |16| as a function of f3. Compared to our previous measurement, the range in 
/3 is considerably increased. Figure |T6| shows that F^^^^ rises as /? decreases. This is expected 
from the QCD evolution of the parton densities in the pomeron. 

The f3 dependence of F^^^^ is sensitive to the dynamics of the 7*-pomeron interaction and 
can distinguish between different pomeron models. In Fig. |1^ ^2^''^^ is compared with the 
predictions of various models. In the model of the pomeron is represented by a 

single gluon leading to photon-gluon fusion followed by subsequent colour compensation. The 
colour compensation is considered to be sufficiently soft so that the dynamical properties of 



the photon-gluon fusion process remain unchanged. The predictions of |Q, shown as the 
dashed-dotted and dotted lines, were normalized to the value of ^2^''^'' at (3 = 0.5, Xp = 0.003 
and = 31 GeV^. The model fails to reproduce the rise of ^2^''^'' towards small (3 values and 
the dependence at large (3. The prediction of the Hard-POMPYT model (full line) where 
the 7*-pomeron interaction results in a quark- ant iquark final state (called the hard component) 
and where the quark momentum density of the pomeron is given by /?/(/?) oc j3 ■ [1 — (3) also 
fails to describe the measured (3 dependence of F2^'^^- However, agreement can be obtained by 
the inclusion of a soft component in the pomeron leading to the form 

/?/(/?) oc/3-(l-/3) + |- (1-/3)^ 
as suggested in the NZ model ||3l]]. A fit to the data yielded g = 0.78 ± 0.32 which is in 



agreement with the NZ prediction of g ^ 1. The fit, shown as the dashed curve in Fig. |1^, 
gives a good description of the data. 



8 Conclusions 

A novel method was used to extract the diffractive cross section in deep-inelastic electron-proton 
scattering. Previous analyses were based on pseudorapidity gap distributions and depended on 
the detailed modelling of the diffractive and nondiffractive contributions. The new method is 
based on the measurement of the mass Mx of the system X resulting from the dissociation of 
the virtual photon and assumes that, for nondiffractive scattering, low lnM|- of the hadronic 
system observed in the detector are exponentially suppressed. The exponential slope and thus 
the nondiffractive contribution were obtained from the data and were found to be independent 
of the specific form of the diffractive contribution. 

The W dependence of the diffractive cross section da'^l^J_^xN{^x,W,Q'^,Mi\f < 4:GeV)/dMx 
measured at large between 10 and 56 GeV^ yielded a value of cT^ = 1.23 ± 0.02{stat) ± 
0.04:{syst) for the pomeron trajectory averaged over t. The same W dependence was found for a 
subset of the events which have Mx < 3 GeV and which are characterized by a large rapidity gap 
and a negligible nondiffractive background. The value of was obtained under the assumption 
that the contribution from longitudinal photons is zero. Assuming af^^ = {Q"^ /Ml)(j'^^^ leads 
to a larger value of = 1.28 ± 0.02(stat). The value for cF^ measured in this experiment is 
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substantially larger than the result found for the soft pomeron in hadron - hadron scattering 
averaged over t, = 1.05, a value which is also consistent with data on p° production by 
real photons at HERA The observation that is substantially larger in diffractive 

DIS than expected for the soft pomeron suggests that in the kinematic region of this analysis 
a substantial part of the diffractive DIS cross section originates from processes which can be 
described by perturbative QCD. 
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Mx 

range 
(GeV) 



range 
(GeV^) 



W 

range 
(GeV) 



obs 



J\fegas 



J^nondiff 



da/dMx istat ±syst 

(nb/GeV) 



< 3 



10-20 



60- 74 
74- 90 
90-110 
110-134 

134-164 
164-200 



49 
67 
87 
91 

97 
118 





24 
8 

8 
24 



2 
1 
1 

4 

1 
1 



± 2 

± 1 

± 1 

± 2 

± 1 

± 1 



36.7 ± 5.7 
46.6 ± 6.4 

55.5 ± 9.2 

68.6 ± 8.4 
78.4 ± 9.0 

94.7 ±11.3 
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2 



±12 
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13.7 
15.9 
20.4 
31.0 
32.0 



± 3.2 
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± 2.5 

± 4.1 

± 4.0 

± 4.4 



3.5 
5.1 
4.0 
2.9 
2.2 
1.8 
6.9 
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2.5 
3.1 
2.6 
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57.5 
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± 6.3 
± 7.4 
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10.3 

8.1 
11.8 
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77 
52 



13 
3 



± 6 
± 2 



23.8 
26.9 



± 3.0 
± 3.4 



± 



1.7 
2.9 
I 1.3 
3.2 



Table 1: Results for the diffractive scattering via 7*p X + N, where N is the proton or 
dissociated nucleonic system with mass Mjv < 4 GeV. The table contains, for each bin, the 
ranges of Mx, Q^, and W; the number of observed events in the bin, N'obs', the number of 
background events from electron gas scattering, A/"^^"* (obtained from the number of events 
found with unpaired electron bunches and scaled up by a factor of ~ 8 to account for the 
difference in electron currents) and from nondiffractive scattering, J\f^°"-'^^ff-^ the value of the 
differential cross section da/dMx at the Q'^ values of 14 and 31 GeV^ and the logarithmic means 
of the W intervals averaged over the specified Mx range and its statistical and systematic errors. 
The overall normalization uncertainty of 3.5% is not included. 
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Q' (GeV^) / Mx(GeV) 


2m,r — 3 


3-7.5 


7.5 - 15 


2m^ — 15 


14 


2.9 ± 0.7 % 


4.8 ± 1.0 % 


4.6 ± 1.2 % 


12.3 ± 1.7 % 


31 


0.9 ± 0.2 % 


2.7 =b 0.6 % 


3.4 ± 0.8 % 


7.0 ± 1.0 % 



Table 2: The ratio r*^^ = j dMx(T%UxNi.MN < 4 GeV)/a*fp integrated over different Mx 
bins atW = 181 GeV and = 14 QeV^ and 31 GeV^, respectively. 
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Figure 1: Illustration of deep inelastic scattering: nondiffractive scattering (top); diffractive 
scattering without breakup of the proton (middle); diffractive scattering with proton dissocia- 
tion (bottom). 
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Figure 2: a) The corrected mass Mx'^'^"^' versus the generated mass Mx for MC events 
generated with POMPYT in the W, range of this measurement, b) The same for the average 
values of M|^'^^'^ and M^'^^'^". The vertical bars show the rms uncertainty of M^*^^"^ for a 
single measurement. 
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Q' = 1 4 GeV' 




5 10 15 20 5 10 15 20 



Mx (GeV) Mx (GeV) 

Figure 3: Left: Mx distributions for MC events on the generator level (histograms) and 
after reconstruction and measurement simulation (points with error bars) as a function of Mx 
for different W bins at = 14 GeV^ for weighted POMPYT (see text). Right: The ratio 
j^MCmeas jj^MCgen measured and generated MC event numbers as a function of Mx- In all 
plots the hashed areas show the Mx regions used for the determination of the diffractive cross 
section. 
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Q' = 31 GeV' 
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Figure 4: Left: Mx distributions for MC events on the generator level (histograms) and 
after reconstruction and measurement simulation (points with error bars) as a function of Mx 
for different W bins at =31 GeV^ for weighted POMPYT (see text). Right: The ratio 
j^MCmeas jj^MCgen measured and generated MC event numbers as a function of Mx- In all 
plots the dashed areas show the Mx regions used for the determination of the diffractive cross 
section. 
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Figure 5: DIS events from data measured in the Q range 10 - 56 GeV^: the scatter plot of 
Mx versus W. The events with rj^ax <l-5 are shown as larger dots; they concentrate at small 
values of Mx- 
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Figure 6: Distributions of Mx at = 14 GeV^: a) - c) for the W intervals 90 - 110, 134 
- 164, 200 - 245 GeV. d) - f) Distributions of lnM|^ for the same W intervals. The shaded 
histograms in a) - f) show the events which have rj^ax < 1-5 corresponding to a rapidity gap in 
the detector larger than 2.2 units, g) - i) Distributions of the rapidity gap Arj for the same 
and W values. The histograms show the distributions for events with Mx < 3 GeV (shaded) 
and Mx = 3-7.5 GeV (skewed hatched). Here Mx is the corrected mass. The distributions are 
not corrected for acceptance effects. 
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Figure 7: Distributions of Mx at Q"^ =14 GeV^: a) - c) for the W intervals 90 - 110, 134 - 
164, 200 - 245 GeV, d) - f) plotted versus In Mjr for the same W intervals as in a) - c). Here 
Mx is the corrected mass. The distributions are not corrected for acceptance effects. Shaded 
histograms show the prediction of the NZ model for diffractive production; dashed histograms 
show the prediction of CDMBGF for nondiffr active production. The MC events were passed 
through the standard ZEUS detector simulation. The CDMBGF distributions were normalized 
to 85% of the data while the NZ distributions were normalized to the observed number of 
diffractive events with Mx > 1.7 GeV; in the NZ model, diffractive events are generated for 
Mx > 1.7 GeV. 31 
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Figure 8: Distributions of InM^ + \n{s/W^) for the W intervals 90 - 110 GeV (dotted), 
134 - 164 GeV (dashed), 200 - 245 GeV (sohd) {laW^ = 9.0 - 9.4, 9.8 -10.2, 10.6 - 11.0) at 

a) =14 GeV^ and b) 31 GeV^. Here Mx is the corrected mass; the distributions are the 
measured ones, not corrected for acceptance effects. For each the three distributions were 
normahzed to the same number of events. 



32 



CDMBGF Monte Carlo 




1 =- 



3 4 5 6 7 

In(Mx^) 




3 4 5 6 7 

In(Mx') 




3 4 5 6 



In(Mx^) 

Figure 9: Comparison of InMx distributions for Monte Carlo events produced with CDM- 
BGF at the generator level (dashed histograms) and at the detector level (full histograms) 
for =14 GeV^. At the generator level, Mx denotes the true mass of the particle system 
produced up to pseudorapidities of 77 =4.3 (forward edge of calorimeter). At the detector level 
Mx denotes the observed and uncorrected mass. The straight dashed and full lines show the 
fits of the exponential slopes to the CDMBGF distributions. 
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Figure 10: Distributions of lnM|- for the W intervals indicated at = 14 GeV^. The data 
are shown with error bars which give the statistical errors. Here Mx is the corrected mass. 
The distributions are not corrected for detector effects. The solid lines show the extrapolation 
of the nondiffractive background as determined by the fits (see text). The dotted histograms 
show the predictions for nondiffractive scattering as calculated from CDMBGF. The CDMBGF 
distributions were normalized to 85% of the number of events in the data. 
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Figure 12: Example of a fit for the determination of the nondiffr active background in the W 
interval 134-164 GeV at =14 GeV^. The data distribution of lnM|- is shown (uncorrected 
for detector effects) with error bars which give the statistical errors. Here Mx is the corrected 
mass. The dotted hne shows the fit performed with D2 — 0. The beginning and the end of 
this line show the InM]^ range over which the fit was performed. The solid line shows the 
nondiffractive background as determined by the fit (see text). 



36 



ZEUS 1993 



> 

c 



,200 



150 



Q' = 1 4 GeV' 
• Mv < 3 GeV 



100 



b 



50 




100 



200 300 
W(GeV) 



> 

(V 

o 
\ 

C 



,200 



150 



Q' = 1 4 GeV' 

• 3 < Mx < 7.5 GeV 



100 



b 

X) 



50 




100 



200 



300 
W(GeV) 



> 

c 



,200 



150 



q2 = 14 GeV' 

• 7.5 < Mx < 1 5 GeV 



100 



b 

XI 



50 



100 



200 




300 
W(GeV) 



50 
40 
30 
20 
10 


I 

50 
40 
30 
20 
10 


300 
W(GeV) 



Q' = 31 GeV' 

• 3 < Mv < 7.5 GeV 




100 



200 



300 
W(GeV) 



Q' = 31 GeV 
^ • 7.5 < Mx < 15 GeV 




100 



200 



300 
W(GeV) 



Figure 13: The differential cross sections da'^''^^^ {j*p XN)/dMx as a function of W 
averaged over the Mx intervals 2171^^ - 3, 3 - 7.5 and 7.5 - 15 GeV at = 14 and 31 GeV^. 
The inner error bars show the statistical errors and the full bars the statistical and systematic 
errors added in quadrature. The overall normalization uncertainty of 3.5% is not included. The 
curves show the results from fitting all cross sections to the form da^^^^ /dMx oc (VF^)(^^jp~^^ 
with a common value of a^, see text. 
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Figure 14: Sensitivity of the value of cF^ to the different sources of systematic uncertainties. 
The central line and the shaded band show for the standard fit the value of and ± 1 s.d. 
The dots give the value with its uncertainty obtained by repeating the analysis for each 
systematic check labeled 1 through 17 as described in the text. 
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Figure 15: The diffractive structure function as a function of Xp from this analysis 

(sohd dots). The error bars show the statistical and systematic errors added in quadrature. 
Also shown are the results from our previous measurement]^] (open dots) and from HI |^ 
(stars) obtained at slightly different /5 and Q'^ values. 
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Figure 16: The diffractive structure function F2 as a function of at x^p = 0.003. The 
error bars show the statistical and systematic errors added in quadrature. The full line shows 
the prediction of Hard-POMPYT. The dashed line shows the prediction of Hard-POMPYT 
with an additional gluon contribution suggested by the NZ model and fitted to the data (see 
text). The dashed-dotted (dotted) line shows the prediction of the pomeron model based on 
the photon gluon fusion dynamics at 14 (31) GeV^ (see text). 
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